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Abstract. In the context of loop quantum cosmology, we consider an inflationary 
era driven by a canonical scalar field and occurring in the semiclassical regime, 
where spacetime is a continuum but quantum gravitational effects are important. 
The spectral amplitude and index of scalar perturbations on an unperturbed de 
Sitter background are computed at lowest order in the slow-roll parameters. The 
scalar spectrum can be blue-tilted and far from scale invariance, and tuning of the 
quantization ambiguities is necessary for agreement with observations. The results 
are extended to a generalized quantization scheme including those proposed in the 
literature. Quantization of the matter field at sub-horizon scales can provide a 
consistency check of such schemes. 
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1. Introduction 

One of the goals of modern theoretical physics is to find a consistent quantum description 
of gravity and the other interactions in Nature. This description should be able to 
reproduce the phenomenology of the standard particle and cosmological models and 
predict new, observable effects. So far, only two candidates seem to fulfil at least part 
of the necessary requirements of such an ambitious program. One is string theory: its 
nonperturbative structure is still under intense scrutiny but some of its contents were 
already exploited to motivate cosmological scenarios where the visible universe is an 
extended four- dimensional object (a brane) living in a higher-dimensional, continuous 
spacetime (see [H [2J, [3J S] for some reviews). 

An alternative, nonperturbative candidate is loop quantum gravity (LQG) [5, 
El d IHJ , where the formulation of the quantum theory is carried out via the 
Hamiltonian formalism. The specialization to isotropic backgrounds, and in particular 
to a Friedmann-Robert son- Walker (FRW) metric, is referred to as loop quantum 
cosmology (LQC) (see [9] for a review). The resulting spacetime is four-dimensional 
by construction. 

The symmetry reduction of the theory is performed at the kinematical level, that is, 
before solving the Hamiltonian constraint [TQl [HI 121 [131 EH] ■ This procedure resembles 
standard minisuperspace quantizations, in the sense that the symplectic structure is 
reduced at the classical level [HI [T5J. However, the identification of the canonical 
variables and their quantization closely follow the formulation of loop quantum gravity. 
The states in the resulting Hilbert space, constructed in [15], are in fact regarded as 
symmetric states of the full theory. 

At very small scales, below a critical scale factor a;, spacetime is discretized and the 
Hamiltonian constraint equation can be written as an evolution equation with discrete 
time [T3l[T6] . Above another critical value a*, larger than the Planck length £p\, classical 
dynamics is recovered, while in the intermediate, semiclassical regime a\ < a <C a* 
spacetime can be treated as a continuum but nonperturbative quantum effects modify 
the standard cosmological evolution. The evolution equation can then be compared 
with the standard Wheeler-DeWitt equation [13l HH HZl HH HH [201 EH [22]. 

Several authors have studied the semiclassical limit of LQC in the context of early- 
universe inflation [111 [IS [IS [23 [1Q [2lS (IE 12^ 

or bouncing/cyclic scenarios [T5l [391 SHI EH S21 H3]. Loop quantum gravity can 
leave its imprint on the large-scale structure of the universe and, in particular, the 
cosmic microwave background (CMB). However, the understanding of inflation in this 
framework is still partial and further study of its testable effects is required. 

In this paper, we compute the main cosmological scalar observables from the 
perturbed background equations using standard techniques (e.g. [B]), assuming that 
inflation takes place only in the semiclassical regime. The scalar amplitude and spectral 
index, as well as the index running, are found at lowest order in the slow-roll (SR) 
approximation, that is, on a de Sitter (dS) background. The scalar spectral index can 
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be blue-tilted and far from the Harrison-Zel'dovich spectrum n s = 1. We consider 
two examples of background homogeneous exact solutions and find that quasi scale 
invariance, as required by observations (see [15] for WMAP experiment), is either lost 
or achieved at the price of a certain amount of tuning of the ambiguities. 

While this work was under completion, other authors have considered the issue of 
scalar perturbations without metric backreaction [38]. We will compare their results 
with ours in due course. 

The paper is organized as follows. In section [2] we briefly recall basic concepts and 
construction of loop quantum cosmology. Several aspects and problems of LQC and the 
present treatment are listed in section [31 The inflationary scalar spectrum and index are 
found in section 0] and calculated for a couple of examples in section [51 The analysis is 
extended to other quantization schemes in section El Section [7| is devoted to discussion 
of these results and future prospects. 

2. Setup 

The basic multiplication operator in homogeneous loop quantum gravity are point 
holonomies, which are S77(2)-elements associated with a single point on the space 
manifold and serve to describe quantum scalar fields. The further requirement of 
isotropy implies that quantum states (in the connection representation) are functions on 
a single copy of 577(2). In particular, the holonomies can be written as hi = exp(/i cri) 
in the spin-1/2 fundamental representation, where I = 1,2,3, T\ = — ici/3 G SU(2), 
are Pauli matrices, and /io = \/3/4 is a quantization ambiguity related to the length 
of the holonomy and fixed to a natural value [Hj (in many papers, /io = 1). The 
connection component c = (VK — ja)/2 is conjugate to the density weighted dreibein 
component \p\ = a 2 , where K is the curvature constant (K = for a flat universe, K = 1 
for a closed one; the Hamiltonian formulation in the open case is treated separately 
[SI SI]); 7 ~ 0.2375 is the Barbero-Immirzi parameter, and derivatives with respect 
to synchronous time are indicated as dots. The spatial volume is then V = \p\ 3 ^ 2 |llj . 
while the symplectic structure is {c, p} = k 2 7/3, where k 2 = 87r£pj and curl brackets are 
Poisson brackets. 

The classical Hamiltonian constraint is 

>/b[ + Wn*t = 0, (i) 

where 7i mat is the Hamiltonian of matter, which for simplicity we assume to be a real 
scalar field: 

ftmat = + a 3 V(<J)) , (2) 

where = a 3 <f) is the momentum canonically conjugate to <ft and V is the field potential. 

Volumes and densities are promoted to bounded operators with discrete spectra 
[T8| HH1 US]. The volume and dreibein operators V and p have zero eigenvalue at the 
origin of the discretized 'time' variable it 6 Z, and their inverses fail to be densely defined 



H 
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c(c - VK) + (1 + 7 2 )f 
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operators. However, the inverse-volume operator V 1 is regular (with zero eigenvalue) 
at t — 0: this way the big bang singularity is healed by purely geometrical effects (see 
also pj). 

The inverse volume can be written in a convenient way by using the classical identity 
a~ 3 = [3(k 2 7/) _1 {c, |p| z }] 3 ^ 2-2 ^ [49], [25] and replacing the connection components with 
holonomies, 

3/(2-2/) 



a' 3 



(3) 



^fi lj(j + l)(2j + 

which does not involve inverse, classically divergent powers of a if < I < 1. 

The ambiguity parameter I determines the initial slope of the effective geometrical 
density. A natural choice, often used in literature, is / = 3/4 [26]. To preserve coordinate 
invariance when quantizing geometrical densities, I must be discrete, Ik = 1 — (2&) -1 , 
k G N [HH EE]. Hence one can select the bound 1/2 < Z < 1, which is also favoured 
phenomenologically [25] . 

The other ambiguity j is a half-integer setting the value of the scale factor at which 
classical cosmology is recovered. It arises for arbitrary spin-j representations (in which 
the above trace is taken) of the holonomy [T8J 09]. The fundamental representation 
corresponds to j = 1/2. 

Promoting Poisson brackets to commutators of operators, one defines an inverse- 
volume operator. If the wavefunction of the Universe does not oscillate at small scales, 
difference operators can be approximated by differential operators which also appear in 
the standard Wheeler-DeWitt equation. The matter Hamiltonian, however, is not as 
in standard quantum cosmology, equation ([2]), because of the promotion of the inverse 
volume to an operator dji = a~ 3 with eigenvalues djj = Dia~ 3 , where 



(q + l) l+2 -\q- l\ l+2 (q + - sgn(g - l)|g - l| m 



-q- 



3/(2-2/) 



1 + 2 l+l 

(4) 

here, q = (a/a*) 2 , a 2 = a 2 j/3 and a ; = ^/t/2o^pi- Equation (jlj) is valid for large j and is 
a good approximation already at j > 10. A large j extends the period of super inflation 
|26j, although this does not rule out values 1 < j < 10 2 . As an upper bound, we 
note that no quantum gravitational effects have been detected so far in accelerators for 
energies below E ~ 10 3 GeV. Therefore the characteristic scale should be (v^pi) -1 < E 
and j < 10 30 . 

djj has a maximum at w a*. One recovers classical behaviour when q ^> 1 (Di — > 1), 
while in the limit q 1 (that is, in a regime where LQG effects are important) [IHl [18] 

D^cla 2 ^, (5) 



3/(2-2/) / a \ 3(I-2)/(I-l) 



J + 1/ \ a * 
In the first line we have defined the constant coefficient c v and 



(6) 
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where 

a = aDl l/2 . (8) 

We will often encounter derivatives with respect to In a. This comes out naturally from 
the 'intrinsic time formalism' of standard quantum gravity (as well as in the continuum 
limit of LQC, where the matter Hamiltonian is quantized differently p3]), where In a is 
the time variable of the Wheeler-DeWitt equation. 
In the LQC regime 
a u 

a^—, (9) 

v ~ - 4 ~ * , (10) 
2(1-/)' K ' 

while in the classical regime a ~ a. Most of the equations below, unless stated otherwise 
(with a symbol ~ as in equations (jUJ) and ffTUl) ). are found for any smooth function Di 
and v. 

It will be useful to note that the classical limit is achieved directly from equation 
([5]) when v — 1 (I = 2, c\ — 1). The LQC regime is valid when v < —2, and the 
eigenvalues of the inverse volume operators are well defined even when v — > —2 (/—►()); 
when I = 3/4, v w —13/2. Hereafter j > 3, so that ai < a* and the semiclassical regime 
is well defined. 

2.1. Equations of motion 

The semiclassical Hamiltonian is [18] 



H = a(a 2 + ^) + ^i^J + «(V0) 2 + a 3 V, (11) 

where now 

K* = djJ<j>, (12) 

and for later convenience we have included the gradient term. The inverse-volume 
operator and gravitational Hamiltonian are constructed through holonomies in the spin- 
j and spin-1/2 representation, respectively. The constraint 7i = is the Friedmann 
equation, while combining the Hamilton equations = {0,?^} and tt^ = {ti^^TL} one 
gets the dynamical equation for the scalar field: 

k 2 ( 4> 2 , \ K 

V2 



H< = ^[^ + V\-^, (13) 



«S-^VV + (l + 2iO*ty + D,V^ = 0, (14) 

where H = a/ a is the Hubble parameter, and V is differentiated with respect to <fi. 
At high energies the effective Klein-Gordon equation is that of a canonical field with 
potential W = j DiV^ d<p in a contracting universe with — (l + 2i/) > 3 dimensions. This 
simple remark eventually warns us that the resulting inflationary observables may be 



Inflationary scalar spectrum in loop quantum cosmology 



6 



at odds with experiments, since the anti-friction term feeds the kinetic energy making 
it dominant over the potential energy [23]. 

The equations of motion can be recast in terms of an effective perfect fluid [40J: 

p eff ^ + ^, (15) 
_l + 2z/0 2 

Pcff = ^2A" r ' (16) 

so that 

H 2 = j PcS , (17) 

p cfr + 3H (Pcff + Pes) = . (18) 

In the limit v — *■ 1 these equations reduce to the Klein-Gordon energy density p and 
pressure p. 

From now on we move to the special case of a flat FRW background, K — 0. 
2.2. Quantization schemes 

As we shall see, ambiguities can in general be observed and fixed experimentally 
[25l l26l [49J. While I and j appear at the kinematic level, other ambiguities arise when 
changing the quantization of the Friedmann equation. It is straightforward to extend 
the above equations of motion accordingly. 

Following the classification of [26], we denote with HAM(n), n > 0, the scheme 
where arbitrary positive powers of a 3 djj are inserted into the Hamiltonian. Classically 
these insertions are equal to unity, but upon quantization they become dynamically 
nontrivial. In particular, in HAM(n) both the matter effective energy density in the 
Friedmann equation and the definition of in terms of the momentum tt^ are multiplied 
by Df: p c ff(7r^,0) — > [n]^, 4>) = D^p^D^^^cf)). The above equations are in the 
Ham(O) scheme. 

A second possibility is the Fried quantization [23], where the squared Hubble 
parameter is promoted to an operator and the Hamiltonian is no longer regarded as 
the primary object. When taking expectation values, the Friedmann equation reads 
H 2 = {H 2 ) oc (djj) (Timat)? and there appears an extra factor D\ relative to the Ham(0) 
Friedmann equation. The scalar momentum and field equation are unchanged. 

A third case (which we will call Pleb) occurs when the classical gravitational 
Hamiltonian is derived from the self-dual Plebanski action instead from the Einstein- 
Hilbert one [36]. Then 7ipi e b = T~t/s/p classically (corresponding to a time 
reparametrization dt — > y/pdt), which upon quantization gives extra factors of D]^ . 

Equations ([TBI and ffl4l) can be generalized to all these schemes as 

D?1+ i H 2 = + Df-+Vj , (19) 

D n 4 +l 



I 



a 2 



-V 2 + [(1 + 2u) + 2n 3 (u - l)]H<j> + £>f" 2+ % = 0, (20) 
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Scheme 


m 


n 2 


% 


b 2 


Ham(w) 


n 


n 


n 


(v-l)n 2 + v 


Fried 


-1 








v-l 2 + u 


Pleb 


1/3 


1/6 


1/3 


(4i/ + 5)/3 



Table 1. Ambiguities in different quantization schemes. The parameters 6i and 62 
are defined in Sec. [5] 



where n\, n 2 , and % are shown in table [TJ Here we have allowed also for a general 
correction to the gradient term, which was introduced in |28J for n 4 = (equation (11 ID ) 
on phenomeno logical grounds. Later on we shall see how to constrain it. 

At first we shall consider in detail the case Ham(O). It is interesting to 
check whether the results thus obtained are qualitatively robust across inequivalent 
quantization procedures, which is the case as argued in [26]. The effect of the 'secondary' 
ambiguities n x is discussed in section [61 



2.3. Superinflation and slow-roll parameters 

Using D1/D1 = 2H{1 — v\ one finds the Raychaudhuri equation: 

K 2 j) 2 

H = -(2 + u)^. (21) 
6 Di 

Since v + 2 < 0, the universe superaccelerates in the semi-classical phase (H > 0). The 
first parameters of the slow-roll tower are chosen as in the classical case: 

(22) 
(23) 

(24) 

The first expression is purely geometrical, does not depend on the effective action, and 
has the advantage to define inflation precisely when e < 1. According to the second 
formula, the condition \r\\ -C 0(|i/|) is sufficient to neglect <fi with respect to the modified 
friction term in the Klein-Gordon equation. 

From equations f[T3"j) and fl2T|) . the evolution equations of the SR parameters are 

v 




2He 



T]-{l-u) + 



(25) 



2H(2 + u)_ 

& 2He[e - r] - (1 - u)], (26) 
r) = H(e V -e)- (27) 
The extra factor (1 — u) in equation ( f26l) will contribute to modify the scalar spectrum. 

Since we want to extract the cosmological observables at lowest order in the SR 
parameters, it is sufficient to assume a de Sitter background (H = const). Note that in 
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the traditional SR approximation all the observables are written at the same SR order, 
but for our purposes the dS computation will be enough, even if it does not allow such 
a systematic truncation (see section l5\3|) . In particular, on can ignore all derivatives of 
H and V in the equations of motion and metric perturbations in the linearly perturbed 
equations. 

3. Caveats and open problems of semiclassical LQC 

Before proceeding, we underline several issues characterizing different formulations and 
approximations of loop quantum cosmology. 

3.1. Conformal rescaling 

In a flat FRW background, the scale factor can be rescaled by a conformal transformation 
which leaves the classical equations of motion invariant. If the spatial homogeneous 
slice is noncompact (K = 0), all integrations must be performed over a fiducial cell, of 
size ~ a, within that slice. Then the fiducial cell also has the above 'gauge' freedom 
(in a closed universe, one can fix the cell to have unit volume). On the other hand, 
the characteristic scale a* is fixed by the theory and the ratio a/a*, as well as its 
functions and the statement that 'the semiclassical regime is such that a < a*', are 
not conformally invariant ([IS], appendix B.2). Then the theory and its observable 
implications would depend on the choice of a, i.e. on specific initial conditions, even 
in the recent '/I quantization' [T5], [52]. As we shall see, although the normalization of 
the scalar spectrum depends on a*, the spectral index does not, and contains only the 
ambiguity v. However, both observables are derived in the semiclassical limit and the 
gauge dependence is, at best, implicit. 

An escape route from the rescaling issue might be to start with a closed universe, for 
which there is no ambiguity in the scale factor. If the background is superaccelerating, 
the classical connection is c = a since the first term rapidly 

vanishes. Then the universe loses memory about the initial curvature and one can use 
the equations in flat FRW. However, one may not be entitled to rely upon this picture 
too close to the singularity, when the scale factor, although expanding, is still very 
small; this might indeed happen in semiclassical LQC inflation. Also, there may be 
additional complications when perturbations are taken into account. This point should 
be addressed in the near future. 

3.2. Improved formulations of LQC 

The previous Hamiltonian constraint equation flTTj), in the spin-1/2 representation and 
in the flat case (K = 0), was generalized through a WKB approximation directly 
at the level of the difference equation [20l ED] and via a path integral quantization 
[291 [32]. Corrections to the classical constraint are generated by the discrete structure 
of spacetime when smoothing the dynamical difference equations near the Planck scale 
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to the continuum limit. Upon quantization of the connection in terms of holonomy 
operators, the factor c 2 in equation ([T]) is modified as c 2 — ► /!q 2 sin 2 (p c) [311 [31] and 
the relation between c and a becomes a = — (7po) _1 sin(/i c) cos(poc). The effective 
Friedmann equation is then 

# 2 oc p cff (l - V (28) 



Pcrit , 

where p cr j t is a critical density. This kind of evolution has applications for bouncing 
scenarios [43] . 

At the time of completion of this paper, we assumed that inflation takes place at 
energies p e fr > Pcrit, that is, when p c <C 1; then the standard Friedmann equation 
(I17p can be used as in the rest of the discussion. A problem might have been that 
recollapse can in principle occur during the inflationary regime [31] . since p crit can be 
made arbitrarily small by increasing the scalar field momentum. However, the above 
'po quantization' does not take into account the expansion of the fiducial cell in a FRW 
background. This was done in the novel p quantization described in [15l [52] , where 

p -> p oc — . (29) 
a 

The functional form of the inverse volume eigenvalues turns out to be modified; the 
Friedmann equation is still equation f[28l but with p crit fixed and around 80% of the 
Planck density. The equation for the scalar is the classical one after neglecting rapidly 
damped quantum corrections ([15J, appendix B). One would get a standard Friedmann 
equation when pc <C 1, that is, roughly when Ht-p\ <S 1- This corresponds to the 
classical regime, and all quantum corrections are neglected accordingly. Such a picture 
is qualitatively different from the semiclassical setup enunciated so faiyi 

There is another, independent source of concern. The parameter j is related to 
representation choices when quantizing inverse volumes in the matter Hamiltonian. A 
similar representation ambiguity arises also in the gravitational sector. Thus the classical 
term y/\p\ becomes a function Sjj which is approximated by a when a > a* (I is fixed 
to 3/4 [321 [M] or 1/2 [IS]). 

In order for the inverse volume corrections to have significant impact on the 
dynamics, one has to quantize the Hamiltonian using the representation of the gauge 
group different than fundamental (j > 1/2). As a matter of fact, the semiclassical regime 
discussed above makes sense only in this 'large j' representation; otherwise, there is no 
intermediate stage between the discrete quantum regime and the continuum classical 
limit. 

For consistency, we should choose the same representation both for geometric 
objects in the gravitational term of the Hamiltonian and for those in the matter sector. 
However, in the higher j case the Hamiltonian constraint is a difference equation of 
higher-than-second order. This may lead to an enlargement of the physical Hilbert 

| Also, in the p, version of LQC the quantum connection c(a, K) is not a linear function of K |53j . The 
role of the curvature in expanding and bouncing scenarios is still to be assessed. 
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space and, as a consequence, to the presence of solutions with incorrect large-volume 
limit [32]. Even if this were not the case, there is evidence (in 2+1 dimensions the proof 
is actually complete) that LQG has a well-defined continuum limit to quantum field 
theory only in the fundamental representation of the gauge group [55]. As said above, 
for small j the semiclassical phase either is reduced or disappears altogether. 

These issues are not completely settled but they would deeply modify the 'old' 
quantization adopted here as working model, with corresponding change in the physical 
picture of the early universe^. Bearing them in mind, we shall try nonetheless to compute 
the power spectrum and the scalar spectral index. It would be interesting to compare 
our results in the absence of metric backreaction with those with full perturbations 
(initiated in [57]; the same background equations are used). This would allow to verify 
what kind of approximations in the treatment of the cosmological problem preserve the 
physical features of the model. 



3. 3. Large j and the number of e-foldings 

One can evaluate the maximum amount of superinflation during the semiclassical regime 
by requiring it to start at the Planck era (a = a(ti) = a{) and to end at a(£*) = a*. The 
number of e-foldings N(t) = ln[a en( j/a(t)] is then 

iV max = ln(a,/ai) = ~ ln(j/3) . (30) 

Therefore, j ~ 10 50 in order to get iV max w 57, which is far beyond the allowed upper 
limit determined by particle physics. Then semiclassical superinflation alone is not 
sufficient to solve the flatness problem. One might invoke some mechanism lowering 
the minimum number of e-folds (an example in classical cosmology is thermal inflation 
[551 [55] ) but j would be still very large; for example, iV max m 25 requires j ~ 10 22 , 
unrealistic in the quantum theory. 

Assuming the universe accelerates even after the semiclassical-to-classical 
transition, would it be possible to see effects of an early semiclassical stage in the 
power spectrum? The range of observable large-scale structures roughly spans the last 
ten e-foldings (e.g., [HD]); these will pertain the classical period, and only predictions of 
standard inflation should be expected. The details of the model may still be tuned so 
that to get a semiclassical modification of the scalar spectrum at low multipoles (very 
large scales); in particular, any deviation from scale invariance in the calculations below 
should be thought of as occurring in that region of the spectrum. 

Moreover, we note that, as argued in [61], the horizon and flatness problems are 
solved if the relative change in the comoving horizon (aif) -1 , not that of a" 1 , is large 
enough. Namely, inflation is more precisely characterized by the number of e-foldings 

§ For recent insights on them, we refer the reader to |56j . There, both the large- j representation 
problem and the matter/gravity consistent representation problem are reinterpreted, and eventually 
relaxed, as natural features related to those of inhomogeneous backgrounds on a rigid lattice. 
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In the standard case H is almost constant and this definition coincides with the previous 
one. In the semiclassical case, it would enhance iV max by a term lii(H*/H{), which is 
positive since H increases in time. Hence, a smaller j is required to get a certain number 
of e-folds. 

In the case under scrutiny, we will compute the power spectrum under the 
assumption of a de Sitter background, later extending the result to a situation where 
the Hubble parameter can vary in time as {H^ 1 )' ~ 0(f) (see section 15.31) . Then 
superacceleration can in principle relax the high-j problem. 



4. Perturbed equations and observables 

We perturb the scalar field around the homogeneous background 4>o(t), and substitute 
0(x, t) = 4>o(t) + <50(x, t) in equation (fl"4l . Going to momentum space, one has 



5k + (1 + 2v)H5^ + 



D,k 2 



0. 



(32) 



where H, Di and V depend only on time and a bold subscript indicates the kth Fourier 
mode. 

The dS approximation (H const) is self-consistent if 

|e|,M«o(M); (33) 

\n\ <C 0(\v\) is required to neglect the mass term in equation (1321) . proportional to 



H 2 



;i + 2z/)(e + r / )-7 ? 2 -e 2 -^ 



(34) 



with respect to the mass given by equation (139|) below, ml oc H 2 v{y+Y). For \v\ — > 0(1), 
these conditions are the usual SR constraints, which are too restrictive in this scenario 
when \v \ is large. 

The slow-roll condition |e| 0(|f |) can be restated as <fi 2 D\V . Thus, in general, 
if the kinetic term dominates over the potential one has |e| ~ 0(|v |) (also, D\ <^ 1 in 
the semiclassical regime). Then one may wonder whether the anti-friction term easily 
causes the kinetic energy to overcome D\V and exit the slow-roll regime. Later on we 
shall discuss what happens if equation (133]) is relaxed to |e|, \r]\ ~ 0(|z/|), and argue that 
e = constant is a sufficient consistency condition. 



4-1. Mukhanov equation 

It is convenient to change time coordinate and define a generalized 'conformal time 



dr.* 
a 



(35) 



|| The name 'conformal' is technically incorrect here, since the metric with this time coordinate is not 
equivalent to the Minkowski metric up to a conformal transformation: ds 2 = —a 2 dr 2 + a 2 dx 1 dx- 1 . We 
shall keep it to distinguish from synchronous time t. 



Inflationary scalar spectrum in loop quantum cosmology 



12 



which reduces to the standard conformal time in the classical limit. In the following, 
primes will denote derivatives with respect to r. In analogy with the general relativistic 
case, one can define a Mukhanov variable «k = a <50k- Then equation ( |32|) . ignoring the 
term in V, become: 



Uu' k + (k 2 + ml) w k = 0, 



where k 2 = Ikl 2 and 



m. 



aH(u 



a 



a: 



= -— + -/, 



a a 
aHv 



H' v' 
aH{v+l) + — + — 
H v 



(36) 

(37) 
(38) 

(39) 



we have not yet assumed that H and v are constant and these expressions are exact. It 
is straightforward to see that in the classical limit (A = const) f v = 0, m 2 = —a" /a, 
and one recovers the standard Mukhanov equation. 
Some useful formulae to get equation (l36il are: 



(40a) 
(406) 
(40 c) 
(40 d) 



d 


1 d 


d 2 


1 /d 2 


a' d 


dt 


a dr ' 


dt 2 


a 2 \ dr 2 


a dr 



a' 

w u 

a 



5<f>'=- 
a 

5<p" = - 
a 

— = ativ , — = — 



u 



2—u' + 



a 



.ot 



,2 



a ' 



2 —- — 

a z a 



u 



a: 



a 



v 



a' H' 
2- + — + - 
a n v 



a a a 

Equation ( 1361) can be solved analytically only for very special functions f v (r) and 
m 2 (r) and at this point we limit the discussion to the LQC de Sitter regime with v = 



const, where one can exactly integrate equation (1351) via equation (jHj) and get 

1 



(41) 



t is an increasing function of synchronous time, running from (t = 0, vanishing scale 
factor) to +oo (t — ► +oo). Then 



nr, ^ - I - + 1 ) — 



(42) 
(43) 



in agreement with equation (30) of [28] . 



% One can also use the standard conformal time dr = At /a and the above definition of to 
get a Mukhanov equation u'^ + (Dik 2 — a" /a)uk = 0; its solution is (a linear combination of) 
\uk\ oc y/—kT Z± ll {c v v~ 1 H* J ~ 1 kT v ) for any Bessel function Z M of order [i = (2v)~ 2 + 1 + v~ x ( 62J , 
formula 8.491.4). Again one obtains the power spectrum found below [38]. Also, the definition of 
horizon crossing D;fc 2 — a" /a « coincides with that given below, since a" /a = v(y + l)(aH) 2 in dS. 
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The solution is 

«k = ^(-^^^(-fcr), (44) 

where B v ^ is a constant and H^ 1 ,, 2l/ -, is the Hankel function of the first kind of order 
1 + l/(2u) (see [62], formula 8.491.3). We shall use the two asymptotic limits 

(£)-", x<<1 , (45 ) 
HP(x) ~ \ — e'("-5"-5) , x > 1 . (46) 

^ V 7TX 

where in the first line \i > is assumed. 

In the long wavelength limit k/(aH) — ► (— for — > 0), when the mode with 
comoving wave number k is outside the horizon, one has 

21+1/(2*) / i \ 

|u k | ~ B Vtk T[l + 1/(2//)] (-—) . (47) 

Here we have use the positivity of 3/4 < 1 + l/(2v) < 1 to drop terms proportional to 

At early times (short wavelengths, kr —>■ +oo), perturbation modes are still inside 
the Hubble horizon and asymptotically we have 

^2 



«k ~ B Vtk J- (-krf-^ 2l/) e- ikT , (48) 

up to a phase dependent on v. When v — 1 (large volumes), one recovers the Bunch- 
Davies vacuum solution (positive- frequencjQ plane waves in the Minkowski limit H — > 0) 
«k ~ e~ lkT / \/2k if, and only if, = y / 7r/(4A;). In general, the friction term in equation 
( f36|) does not vanish and the evolution of the particle perturbations is not driven by 
a wave equation. This is expected in the LQC framework, since small scales are those 
mainly sensitive to loop quantum effects^- Nevertheless, we show below that in the 
quantization scheme Ham(1/2) the Bunch-Davies vacuum arises naturally. 

We can determine B vk by imposing standard commutation relations for the 
quantized scalar field. These read 

0(xi, t), 0(x 2 , t)] =0, (49a) 

[tt (xi, t), 7r (x 2 , t)] = 0, (496) 

(xi, t), 7r (x 2 , t)] = ic5 (3) ( Xl - x 2 ) , (49c) 

where hats denote operators and the conjugate momentum 7r^ is given by equation (|T2"|) . 
In general, the composite operator fr^ = a 3 is not isomorphic to the composition of 
operators djf(f). However, in the semi classical regime the states of the Hilbert space 

+ This justifies a posteriori the choice made in equation to set equal to zero the coefficient of the 

(2) 

other independent solution -ffi +1 /( 2 „) ( — kr). This coefficient may be actually proportional to (1 — v), 



but we will keep equation (|44[) for simplicity. 
* There are other theories where the adiabatic vacuum is not the most natural one, 
noncommutative and trans-Planckian models [63] . 
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on which the operators are defined can be decomposed into a gravitational and matter 
sector, |\&) = a ls rav )a ® |mat)^. Then geometrical and matter operators commute, 
and we can take the eigenvalue of the formers in equation (149 eft . In other words, 
gravitational variables are treated as background while matter fields are quantized^. 
In momentum space and time r, the nonvanishing commutator for the variable «k 

is 

=i^r 1/2 ^ (3) (ki-k 2 ). (50) 

The Mukhanov variable is expanded in terms of creation and annihilation operators, 

u k = w k a k + w* k al, (51) 

where * indicates complex conjugate, w k is the classical solution previously denoted by 
Mk, and 

K 1 ,a| C2 ]=5( 3 )(k 1 -k 2 ), (52a) 
[flki, ak 2 ] = = [4 i; ajj . (526) 
Plugging equation (15TI) into (!50|) and using equation (152 aft , one gets w^w^ — w* k w' k = 

^ i<2 _____ 

i-D z . The small-scale solution ()48l) satisfies this constraint only if 



2kD l 



This is indeed a constant in r, as one can check from equations ([5]), Qj and fj4T]jf f 
The spectrum of the scalar field fluctuations is 

^^<i*i 2 >L,- < 54 > 

where angular brackets denote the vacuum expectation value of the perturbation on the 
state defined as <2k|0) = Vk. Conventionally, perturbative modes cross the Hubble 
horizon when k 2 = k\ pa —m 2 , = v[y + l)(aH) 2 from equations (1411) and (14"5|) . The 
positive coefficient in front of aH can be included in the overall normalization of the 
spectrum, and we can set the wave number at crossing as 

K = aH. (55) 

The solution plugged into equation (154"]) is that at long wavelengths, equations fT4T|) and 
(1531) . although it is evaluated at horizon crossing. Then, up to a numerical factor (equal 
to 1 when v = 1), 

n = 4^(f)\ (56) 



In general, geometrical and matter operators do not act separately on physical states because solutions 
to the Hamilton constraint already incorporate correlations between the two sectors [64 . So operators 
on such states are in general complicated, entangled observables. The issue is delicate and will not be 
dealt with here. 

ff This choice differs from that in ,28] ; see equation (45) in that paper. There, the different normalization 
together with the ansatz fc* = aH for horizon crossing lead to V^, oc H 2 . 
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where Di is given by equation (jSJ) with constant v. Clearly, this quantity is not constant 
in dS due to the presence of D[. 

The ansatz ( 1531) will be crucial to get the correct expression of the spectral 
index. Horizon crossing occurs at scales much smaller than in the standard scenario, 

— 1/2 ~ 

k* ~ D t (aH) ^> k = aH. We note also that, while in the standard case comoving 
modes appear to exit the horizon, (d|/c| _1 /dt < 0), in this scenario they enter it, since 
the comoving effective horizon expands, d|A;*| _1 /dr > 0. This raises an important 

interpretational issue [3H], the limit kr — > being not a natural consequence of the 
inflationary evolution. However, this does not affect the power-law solution below, which 
we shall take as a viable LQC example. This solution is not de Sitter (H ^ const) but we 
will argue in section 15.31 that it can be taken as a background for the dS perturbations, 
as happens in standard inflation. Bearing in mind that the horizon-crossing problem 
still remains, for the moment we let it aside and push further the analytic inspection 
of the model. A posteriori this allows us to find one puzzling result which might be 
relevant in future developments of quantum gravity: namely, quantization schemes may 
be constrained by requiring consistent (pseudo) plane- wave perturbations at small scales. 



4-2. Comoving curvature perturbation 

The spectrum of scalar perturbations is the two-point correlation function of the 
curvature perturbation on comoving flat slices, defined as in classical gravity 

ft=fty, (57) 

to linear order. We are going to show this result by aid of the covariant approach, 
first introduced in [651 ES [67] and recently refined in [681 EH]. The starting point is 
a inhomogeneous spacetime with metric g^ v (signature ( — h++), Greek indices from 
to 3) and filled with a perfect fluid characterized by the energy-momentum tensor 
Tfi V = (p + p)u^u u + pg^u, where = dx^/dt is the fluid four-velocity (m m u m = — 1, 
= 0) and t is proper time along flow lines. For a Klein-Gordon scalar field with 
p = — V (7 0V fJ 0/2 + V and p = —'V a <j)V 7 <f)/2 — V, the energy- momentum tensor is 
Tfj, u = V ^Vvcj) + P9fj.u, which compared with the above expression gives 

u V ^ == _YA m 

which is orthogonal to <fi = const hypersurfaces. Here, ' — d t — u a ^7 a . One introduces 
iV = - / dtV a u a (59) 



dt H, (60) 

where V CT is the covariant derivative and in the last line we have specialized to a 
FRW background; on it, N(t) is just the number of e-foldings. However, here N(x,t) 
and H(pc,t) are time and space dependent. The variable of interest is the covector 
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11^ = —D^N, where = (S^ + u^u u )V u = V M + u^d t is the spatial gradient. Via 
equation fl58l) . one has 

= -d^N + -fd^ , (61) 

where all variables are functions of both time and space. This quantity is the nonlinear 
generalization of the perhaps more familiar linear comoving curvature perturbation 7Z, 
and satisfies a conservation equation at all scales thanks to the conservation of the 
energy momentum tensor, V^T^ = 0, without use of the Einstein equations. At linear 
order, 1Z\ ~ d{JZ, where the index i runs over spatial coordinates; in the flat gauge 
5N = 0, 1Z is given by equation (1571) . with H and now coinciding with the background 
homogeneous quantities. 

In loop quantum cosmology, a covariantly conserved effective stress-energy tensor 
for a perfect fluid is 

T ^u = iPeS + Peff) UpUv + V&Qpu , (62) 

with peg and p e s given by equations ( TT5l) and (TIE]) in FRW for the scalar field. Then one 
can start from a classical Einstein-Hilbert action and get the Friedmann equation ffTTI) 
as the 00 component of the Einstein equations G^ u = k 2 T^. Equation ( |62l) corresponds 
to 

Tf v = |^ V^V,0 + p eS9lxu . (63) 

Equation (163]) is covariant if D\ is regarded as a function of the volume V. The simplest 
way to implement the covariant formalism in LQC is to promote all the above quantities 
(including D\) to time- and space-dependent variables in a locally defined 'separate' 
universe [70]. It is not obvious and likely nontrivial to show that this construction is 
compatible with loop quantum gravity; we shall assume it as working hypothesis, leaving 
its check for the future. In this case, we get again equations f]58]) . flST]) and fl57j) . 

A remark is in order. Within the same framework, one can introduce the uniform 
density perturbation (at linear order) ( c g defined on uniform density slices where 
<5p e ff = = SpcQ. This perturbation will be different from the uniform density 
perturbation ( on slices where 5p = = 5p, and while the latter is equivalent to 
1Z (up to a sign and at large scales), the former is not. The reason is that, for a 
canonical Klein-Gordon field, 5p — 5p = 2V^5(f). The left-hand side vanishes on uniform 
density slices, implying that these are also comoving (5(f) = 0). On the other hand, 
8pcff — 8peE = 2V<j,5(j) + 2(1 — v)(f)(5(f) + C(p)/3, where the term in C comes from the 
perturbed metric in the kinetic term g^d^^d^^. If v ^ 1, in general the right-hand side 
is nonvanishing on comoving slices. Then, TZ fa —( and ( e g are inequivalent candidates 
for the observable spectrum. We skip further subtleties inherent to this choice (including 
the question about conservation of these quantities) and opt for the comoving curvature 
perturbation. The author of [28] chose the density contrast 5p/ p. 
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4-3. Scalar spectrum and index 
The scalar spectrum is 



?k 3 

4 s !?"*'" 



k k 



k| 

2 

^ • (646) 



2H 

5j 

We have chosen a normalization which allows a transparent comparison with the 
standard result; this can be achieved by compensating the normalization of the 
solution for the Mukhanov equation (fixed by commutation relations) with a suitable, 
approximate definition of horizon crossing. Changing the latter would lead to a in- 
dependent factor in front of the amplitude (16461) . which can in principle be tested 
by CMB observations. This would enlarge the parameter space and loosen individual 
constraints on other variables. 
The final result reads 

1 / 1 H 2 ^ 2 



A ° = wA^J ] (65a) 



(2 + u)k 2 H 2 
~ 6(5vr) 2 D f 2 e 

The spectral index is defined as 

din A 2 



n K — 1 = 



d In k 



(65b) 



(66) 



where at horizon crossing din A; = (a/a)dt = uHdt on dS. From equations (165 fej) and 

™ s -l = ^ cl) -l-(l-^)], (67) 

where n$ — 1 = 2r] — 4e is the standard scalar index. Heuristically, since v < the 
scalar index can be blue tilted if the classical one is red tilted. 

This finding agrees qualitatively with that achieved by the 'direct method' of |28j . 
although comparison of the expressions for the scalar index (equation (43) of [2S]) is 
not straightforward. A blue tilt corresponds to a loss of power at large scales. This 
can be understood by the introduction of a cutoff scale, of order £p\, below which no 
modes are produced from vacuum (this is nothing but the 'curvature cutoff' provided 
by quantization of inverse powers of the scale factor). Such a cutoff is ideally present 
also in the standard picture, but here it is rather close to the energies involved (a > £p\). 
A conclusion similar to that of [28J was reached in another minisuperspace cosmological 
model [71] . where the canonical variables and quantum Hilbert space and operators are 
not that of LQG (see [72] for details). Note that also this model has an ambiguity 
analogous to /, which was fixed in [Til 172] . 

However, we will show that the parameter space and freedom in the ambiguities 
allow for both red-tilted and blue-tilted spectral indices, while the scalar index found in 
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[28] is always blue tilted regardless the value of v. Finally, the running of the spectral 
index is 

«s-^ = ^[5e,-4e 2 -e 2 + 4e(l-,)]. (68) 
The last term is 0(e) and may dominate over the others. 



5. Background solutions and the SR approximation 

5.1. Power-law solution 

A power-law solution of the background equations was found in [30] when D\ is given 



by equation (j5J): 

a(t) =tV[(2+")A] j ( 69fl ) 

0(t) =(M (l-)/[(2+.)A] ; (69fc) 

V(0)=K0~ 2(2+,/,V(1 ~ l0 , (69c) 



where (after time reversal) the constant A > 1 in order to have a classically stable 
expanding solution and <p v and V v are appropriate normalization constants. The 
quadratic potential V oc 2 considered in [21] is stable in the LQC regime (A > 1 



for v < —2). The SR parameters read 

e = (2 + i/) A, (70a) 

7] = (2 + u)X-(l-u), (70b) 

e = (2 + z/)A[(2 + z/)A- (1 -v)\. (70c) 

The observable spectral index is constant (a s = 0): 

n s - 1 = 3 ( 1 - - ) - 2 (l + - ) A, (71) 



and is shown in figure[Uas a function of v and A. Most of the parameter space generates a 
strongly scale-dependent spectrum. To highlight a region compatible with observations 
(figure El upper strip) we require —0.1 < n s — 1 < 0, which has the best-fit n s ~ 0.95 of 
WMAP3 [15] as central value. The scalar spectrum can be blue tilted in the semiclassical 
LQC regime but one can tune the parameters to get a red-tilted, almost scale-invariant 
spectrum. In the figures we have considered the conservative interval < / < 1 instead 
of 1/2 < / < 1 (v < -7/2). 

Note that the k~ l actually decreases in time as r, and modes exit the effective 
horizon. 



5.2. Exponential solution 

Another exact solution is of exponential type and can be found either by direct 
calculation or via the mapping found in [30] from the standard solution: 



a(t) = exp ( J^-f ) , (72a) 
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Figure 1. The scalar spectral index (171]) as a function of the parameter A > 1 and 
the ambiguity v. 



A 




Figure 2. Region in the \-v plane predicting a spectral index 0.9 < n s < 1.0 for 
the Ham(O) (upper, dark strip) and Pleb (lower, light strip) quantization scheme (see 
section 16.11) . 



3(l-i/) 



0(t) =0,^ 2 - 1 exp 

2 + v 

v{t) = v v t 2 ^- l \i + A u r s ) ) 

where p and s / 1 are constants and 

2 



- P t S 



3ps 
2 + v 



3 6ps(l - s)c\ 
~k? (2 + z/) 2 

^2 



A. 



2W 



(726) 
(72c) 

(72 d) 
(72 e) 
(72/) 



Reality of the scalar field requires ps(l — s) > 0; setting v — 1 we recover the classical 
solution 
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The SR parameters read 

e (2 + u)(l-s) 
3ps 



t 



(73 a) 



V 




(736) 



2 + u 
3ps 



) 



3ps(s- 1)(1 -v) 
2 + u 



t 



—s 



(73c) 



For positive s, these quantities decay in time and the spectral index is asymptotically 
n s — 1 ~ 3(1 — z^ 1 ); hence 3 < n s — 1 < 4.5, which is strongly blue tilted and excluded 
by observations. 

5. 3. Validity of the de Sitter solution and Lidsey mapping 

The solution equation (jUj) was derived under the assumption that both e and 7] are 
small enough, <C 0(|z/|) (massless field approximation). The last examples, however, 
show that |e| and \r]\ can be 0(|z/|). There may be doubts about the self-consistency of 
the spectral index thus found. 

In standard inflation, e 1 is assumed for the calculation of the dS spectrum into 
which, nonetheless, one inserts a solution (power law) which has e constant of arbitrary 
magnitude (if the constant is small, one has scale invariance). The observables for 
power-law inflation are not at odds with the dS calculation, since the latter agrees with 
the lowest-order perturbative calculation carried out under the assumption e, rj — const. 

As in the standard case, an exact solution of the LQC Mukhanov equation 
at next-to-leading SR order can be achieved when e,r] = const, no matter their 
magnitude. By virtue of equation ( 134"1) . corrected with an extra geometric term in 
e from the backreaction of the metric (which is presently out of the discussion), the 
friction and mass terms (l4"2l and (l4"3"j) would be replaced by [2%2(e, T), v) — l]/r and 
\xk{ e -, r l-, v ) ~ xii.tiVi 1 ')]/ 1 ' 2 respectively, where \2 and %3 are constants built with the 
slow-roll parameters. The solution would be a linear combination of (—kT) X2 Z X3 (—kr), 
where Z X3 is a Bessel function, and will have different asymptotic behaviours relative 
to equation (1441) . The spectrum at horizon crossing, to lowest SR order, is unchanged 
anyway (see the following section, equation (155!) with xo — !)• This might be also 
required by robustness of black-hole results, as hinted in the concluding section. 

The key point is that the adopted 'slow-roll' approximation is actually a 'not-too- 
fast-roll' one, defined in equation (1331) . Roughly speaking, the requirement for the dS 
approximation to hold is much looser in LQC than in classical general relativity, in the 
sense that it is consistent even when the field velocity is not very small. This further 
legitimates the choice of solutions which are not de Sitter, even if the calculations assume 
H = constant. 

The SR parameters can be derived also by using the duality of [30], which relates 
solutions {a(t) , H (t) , (j)(t)} of the LQC equations to solutions of standard general 
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relativity with scale factor b(t), Hubble parameter (3 = b/b and a Klein-Gordon scalar 
ip(t), where b = a 3 ^ 2+u \ = (2 + u)H/3, and ip = -(2 + u)<j>/(3y/Dl). The SR 
parameters calculated for the standard solution {b(t), (3(t), i/j(t)} are 

$ 3 



2 + v 
3 



v 



[V + {!-»)]■ 



(74) 
(75) 



v- 



Hence, a dS solution in LQC corresponds to a standard solution with e = and 77 
this partly explains deviation from scale invariance in the examples. 

Also the scalar spectrum fl56|) can be achieved via the above mapping. One starts 
from the dual system of a Klein-Gordon field in standard general relativity. The two- 
point correlation function of the fluctuation 5ip is evaluated at the standard horizon 



crossing, k = aH, and is given by the Gibbons-Hawking temperature, 
H 2 . One argues that 



/? 2 /(2vr) 



oc 



in agreement with equation ( 1561) 



-fc) d A<i^r) oc 



(76) 



6. Cosmology in a generalized quantization scheme 

Having discussed the implications of LQC for the scalar spectrum in the simple Ham(O) 
quantization, we can easily extend the analysis to the other schemes. The background 
and perturbed dS dynamics can be written starting from the general equations ( 119]) and 
(BP]) . In order to preserve the standard continuity and Friedmann equations, equations 
( TT5|) and (fl6l) become 

+ D 2n2 ~ ni V , 



Pcff - 2£) „ 1+a 



Peflf = 



;(6i+6 2 )-i 



2D? 1 



n,+l 



6j ] D 2m-my 



(77) 
(78) 



where 



6 1 = (l- I /)(ni-2n 2 ), (79) 
b 2 = (i/ + 2) + 2(n 2 -n 3 )(l-i/). (80) 
In the usual scheme Ham(0), 61 = and 5 2 = z/ + 2 (see table [p. For HAM(n 7^ 0), 
61, fo 2 < 0, while in Pleb b\ = and 6 2 < —1. The Raychaudhuri equation is 

b 1 H 2 + H = -b 2 -^ TT . (81) 

The evolution equation for e now reads 



e = h + 2H(e - b± 



[l + n 1 ){l-v) + 



2Hb^ 



2H{e-b 1 )[e-r ] -{l + n 1 ){l-v)l 



(82) 
(83) 
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6.1. Mukhanov equation and scalar field spectrum 



In analogy with what has been done so far, one might define a 'conformal' time and 



canonical variable u = a5<p = aD l 



-(i+m)/2 , 



>. However, the physical result does not 



depend on the choice of time, and we can keep the old definition of r without loss of 
generality. 

The gradient term in the Mukhanov- type equation has an extra factor D™ 4 , which 
can be written as a function of r by using the formula Di m cj v [— za£/t) 2 ( 1-1 /^, from 
equations (JSJ) and (iilj) . The friction and mass terms are 

/„ =a£T(i/-l)(l-2n 3 ) 
l-2n 3 



i-1 

v 



m. 



- aHv < aH[2u + (1 - i/)(l - 2n 3 )] 



if' i/ 



1 



1(1- 2n 3 ) 



The Mukhanov equation and its solution are 



where 



Xo 



2X2-1 , 



r 



\,,-i\2 _i_ xf XoXl 



«k = 



l-- n 4 + l->l 
v 



k 



X1 = [c^^Hf-^r-^k } 

Xo 

1 - 2n 3 1 

X2 = = + n 3^ 7T' 

Iv zv 

X3 = > 1 + 7T ' 

and the limits are taken as n\ — > for all i. Horizon crossing is defined as 

Xi* ~ (aH) X0 fc* « a Xo if . 

At large and small scales, respectively, 

r(x 3 )« 



«k ~ B u>n . >k 



>X3, 



"Xi. 



X2-X3 r X2-X0X3 



_^ l ^X2-l/2 r X2-Xo/2 e -iXi-r xo 



(84a) 
(846) 



(84 c) 
(84 d) 

(85) 
(86) 

(87a) 
(876) 
(87c) 
(87d) 



(89) 
(90) 



The momentum conjugate to is = a 3 D ( 1+ni )0 and the commutation relation (1501) 
becomes 

flUW] = iA ni_1/ ^ (3) (ki - k 2 ). (91) 
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Imposing equation (I9~T|) . the integration constant of the small-scale solution is 



-XO) 1-2X2 £)«i-l/2 

k ~2^~^ 2 XlXo 



R 2 * {-XlT* j -t (Ql? s 



which is independent of r if the consistency condition Xo(l - 2x 2 ) + (2ni — 1)(1 — 1/^) = 
is satisfied: 

2(m-n 3 ) + U- (l-2n 3 )n 4 = 0. (93) 

This constrains the unknown parameter n 4 in almost any quantization scheme. 

• HAM(n): It must be n 4 = for any n ^ 1/2, while for n = 1/2 any value of n 4 is 
allowed. 

• FRIED: Only 2/3 < n 4 = 2u/ [y— 1) < 1 is allowed. Since n 4 should be independent 
of the ambiguity u, it is unlikely that the Fried scheme would lead to a consistent 
quantization at small scales. 

• Pleb: Only n 4 = is allowed. 

These results can be interpreted in different ways: 

• The commutation relation (|9"T1) is indeed a selection rule on the possible quantization 
schemes. 

• The dS approximation is not reliable and an expansion in constant SR parameters 
for the full perturbed equations (i.e. including metric backreaction), or a different 
approach altogether, is required. 

• Any allowed quantization scheme in the symmetry reduced theory is legitimate and 
one has to give up either the commutation relation (149 cp or the oscillator expansion 
(ED. 



Due to the evident limitations of our method, we do not try to solve this issue here. It 
will be important to address it when further developments are available. 

In the Ham(1/2) scheme with n 4 = 0, 2\ 2 — 1 = and the short wavelength solution 
is a plane wave, ~ e~ lkr , thus recovering the Bunch-Davies vacuum. This would be 
also true in the general case 2n 3 — 1 = n 4 , after a time redefinition t Xo — ► r, but such 
condition is not compatible with the above selection rules. 

The scalar spectrum Va, is 



^ _ r> ni-l/2 (fe/a) 2 r 2 (x3) o2>a -l 

x ^-(X2+X3) r Xo(l-2X2-X3)+X2 

(H_V , X2+ 2 HX2 - 



k=k, 



(94) 



oc l^-) (a^H^- 1 )^- 1 , (95) 

in agreement with the previous result V^ik) oc D l l,2 {k/ a) 2 {kr)- 1 / v and equation (156|) . 
The spectrum in Ham(1/2) with n 4 = is precisely the standard one. 
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The comoving curvature perturbation was discussed in the Ham(O) case. The 
only modification is in the form of the effective stress-energy tensor (1621) . which is now 
constituted by the energy density (1771) and pressure (1781) with (p 2 — > — V o -0V°'0. 

6.2. Background solutions, generalized spectrum and mapping 

The exponential solution has an unacceptably blue-tilted spectrum and we shall ignore 
it. The LQC power- law solution is 

a =t 1/£ , (96a) 
<P oc £(i-")("i+i)A ? (966) 

yoct 2(l-,)(n 1 -2„ 2 )/ £ -2 j ( Q6c ) 

v = e -(l-i/)(Tii + l), (9Qd) 
e 2 = e[e-(l-i/)(ni + l)]. (96 e) 

When 6i = as in Ham(O) and Pleb, one can extend the duality of [30] as follows: 

8 = jH , (97) 

and an appropriate redefinition of the potential V{4>) — > W{ip) which is not important 
for the discussion. Then the Raychaudhuri equation (IHTT) can be written as the standard 
one: 

$ = - yV' 2 , (99) 

with a Friedmann equation (3 2 cx = ip 2 /2 + W(ijj). This determines e = 62 A in terms 
of the parameter A from the standard power-law solution a ~ t 1 ^ 3 ^. 

It is instructive to see how different quantization schemes affect the scalar spectrum. 
When 61 = 0, then Xo — 1 an d the spectral index reads 

n s -l = i[-26 2 A - (2ni + 3)(1 - u)\ . (100) 

The Ham(O) and Pleb cases are shown in figure [2j The two schemes have similar 
qualitative features. 

When b\ 7^ 0, the relation between the Hubble parameters H and f3 is not algebraic 
and, in order to preserve equations (1951) and (l9"9~l) . one has to define 



fo 2 /r V , r rr2 N & 2 rV A 6l 



/3 = ^ + 6 1 ^) = ^M - ^1 . (101) 

An alternative might have been to modify equation (198]) so that to simplify the mapping 
H 1— >• f3, but that is the only definition allowing to interpret the spectrum (1951) in terms 
of the dual, standard spectrum for a Klein-Gordon field. By the same procedure of 
equation (176]) . one has that oc D[ 11 ^ 1 ^ 2 (3 2 . For all backgrounds with constant e, (3 is 
proportional to H. 
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Since the mapping is between homogeneous backgrounds, the next step is to get rid 
of those terms in equation fl95l) which bear the explicit imprint of the gradient term. In 
the limit xo ~ * 1> the spectrum agrees with that found via the above line of reasoning. 

The full calculation of the scalar spectrum, however, takes into account at least 
part of the inhomogeneities due to quantum fluctuations, and its expression is more 
complicated than what would be naively expected from the sole mapping argument. 
Hence the relation, as formulated here, between standard and LQC scenarios involves 
only background solutions and not observables. 

In the power-law case, (3 = (3At) _1 , H = (et) _1 and equation (jlOip is quadratic in 
e. It yields 



which is real if A > 4&i/&2- From tabled], in the two cases of interest Ham(h 7^ 0) and 
Fried, the solution is well-defined when A > 4n and A > 4, respectively, thus restricting 
the stability condition A > 1. 

At present we do not know how to fix xo m the Ham(1/2) quantization scheme, but 
obviously it would further change the allowed parameter space. However, the tuning on 
the parameters would not be relaxed but shifted to other regions. 

7. Conclusions 

To summarize, we have analysed the spectrum of cosmological observables generated 
during a period of superacceleration in the semiclassical regime of loop quantum 
cosmology, formulated in the 'old' quantization schemes HAM(n), Fried and Pleb. 
At the level of the background LQC equations, the following assumptions were done 
and discussed: flat metric, classical (or small) connection, large j representation of 
geometrical objects in the matter Hamiltonian, different representation of geometrical 
objects in the gravitational part. By perturbing the background equations without 
taking into account the backreaction of the metric (i.e., on a pure de Sitter spacetime), 
the scalar spectrum of power- law inflation results scale-dependent, unless one tunes the 
parameters of the theory (including ambiguities). The short-wavelength behaviour of 
the perturbed solutions can constrain the ambiguities within each quantization scheme. 

This does not mean that loop quantum cosmology is ruled out by observations. 
The first thing one should check is whether the semi-analytic results and the dS 
approximation agree with a full numerical solution of the Mukhanov equation (I36I) with 
Di given by equation (j4]). Another key point is that we have assumed that inflation 
occurs only in the semiclassical regime. In [2U |26] the LQC semiclassical regime (with 
I = 3/4) was regarded as a selecting mechanism for initial conditions before classical 
inflation, with an intermediate, classical, SR-violating phase in between. This transient 
period may be then responsible for loss of power at large scales, and an indirect signature 
of loop quantum gravity if the field reaches maximum displacement from the potential 




(102) 
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minimum at least 60 e-folds before the end of inflation (here we have seen that low- 
momenta suppression of the spectrum can be also a direct LQC effect). It would be 
interesting to extend the analyses of [24, 26J to other inflaton potentials with general 
ambiguity parameter / and j relaxed to even small values. 

On the other hand, if inflation takes place when the discrete structure of spacetime 
cannot yet be approximated by a continuum, one might be able to put nontrivial 
constraints on the theory by considering the full LQC regime and its difference equations. 
However, one would require explicit expressions for quantum observables which have not 
yet been determined. The task seems difficult but worth being explored. 

Another trend of research might try to extend the matter Hamiltonian to a more 
general content. In [37], for instance, the scalar field is assumed to be nonminimally 
coupled to gravity. Scalars not satisfying a Klein-Gordon equation at the classical level 
can produce qualitatively different, and possibly more viable, stages of inflation. We 
note, however, that certain alternative choices presently in fashion, like the Dirac-Born- 
Infeld tachyon, are motivated mainly by string theory, and their application to LQC [75] 
should be regarded as purely phenomenological. 

One of the most remarkable results of LQG is the proof of the Bekenstein-Hawking 
area law by counting the nonperturbative quantum states of a nonextremal black hole, 
a computation which fixes the Barbero-Immirzi parameter [7_6j [77J [TH [T9J, [80j [8TJ [82] . 
In general, it was shown that the area law is valid for all isolated horizons, including 
the cosmological ones, in the limit of large area (-f^ -1 > 10 3 £pi or so). Then the entropy 
of a dS horizon with radius H~ l is S oc H~ 2 . By interpreting the Friedmann equation 
H 2 oc p c fj as the (integrated) first law of dS thermodynamics — dE = TdS, where 
dE oc H~ 3 dp e g is an infinitesimal energy exchange through the horizon [83J EH [85] , the 
dS temperature is indeed the Gibbons-Hawking temperature Th = H/(2ir) [86]. The 
spectrum of a Klein-Gordon field is = T^. 

It will be important to verify the thermodynamical interpretation of the Friedmann 
equation also for small volumes as in the semiclassical regime. In that case, logarithmic 
corrections to the entropy S ~ H~ 2 + log H + . . . should be taken into account. In fact, 
during semiclassical inflation one cannot use the large- volume expression S oc H~ 2 , 
because H^ 1 < a < a, ~ fpi for j < 10 3 (for larger j there may be an allowed 
window but at this stage the interpretation of V$ as a modified Hawking temperature is 
not clear). Considerations on the canonical partition function of a quantum system in 
equilibrium suggest that the above holographic principle might hold in general, also in 
this regime, because the bulk Hamiltonian (and not the boundary one) is a constraint 
for the spectrum of physical states [87] . 




Although the de Sitter approximation (that is, as a test field with no backreaction 
from the metric) is generally a good one in standard inflation, it is not obvious that 
it is as well justified in LQC. Actually, the symmetry reduction of the full theory is 
performed under the very special assumption of homogeneity and isotropy. But even 
small inhomogeneous perturbations can have drastic effects [TT], and quantization of 
an anisotropic Hamiltonian may lead to a rather different set of equations and spectra. 
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Important open issues to be considered in much greater detail are the choice of the 
(vacuum?) state on which to evaluate the two-point correlation function of the scalar 
perturbation, the choice of gauge (is the flat gauge SN = consistent?) and gauge- 
invariant observables, the adiabaticity of perturbations, and the behaviour of the tensor 
spectrum, yet to be determined altogether. The impact of geometrical backreaction 
should be assessed and what has been shown in this paper be checked by a more 
rigorous calculation in order to be considered reliable. If this were the case, the standard 
technique we adopted would prove to be a useful, relatively simple tool of analysis for 
loop quantum cosmology. 
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